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ABSTRACT 


First results on applying CST Microwave Studio in 
combination with a human voxel model and a 
dedicatedly developed equivalent cardiac cycle source 
are reported. Based on first coarse-meshing 
simulation results it may be concluded that it is 
feasible to model the heart in a three-dimensional 
anthropomorphic voxel model of a human being by 
an equivalent fixed dipole volume source. Through 
changing the amplitude and orientation over the 
cardiac cycle, ECG signals may be obtained from the 
chest area of the voxel model. 
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Il.INTRODUCTION 


* How to determine Electro Cardio Graphic (ECG) 
probe positions on the body? 


* How to combine this with radio? 





Need for model 
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I. INTRODUCTION 
Forward and Inverse Problem of Cardiography 


* In electrocardiography (ECG), the condition of the 
heart is monitored through measuring the body 
surface potential; 


e Measurements are typically taken on the torso; 


e Through a model, the body surface potential is 
related to the epicardial potential = inverse 
problem; 





* For solving the inverse problem the forward problem needs to modeled 
= relating epicardial to body surface potential; 


e From the body surface potential distribution, the body surface charge 
and current distributions may be determined. 
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I. INTRODUCTION 


The Forward Problem 


* Surface modeling: Only boundaries between dielectric regions are modeled, 
regions are considered to be isotropic; 


Volume modeling: The only way to incorporate regions of anisotropic 


conductivity (organs, arteries, bones). 








Numerical computation techniques: 


kV/m 





a) suspended in free space b) standing on pertect ground c) current density 





Finite Elements (FE)* Finite Difference Methods (FDM)** Finite Difference Time Domain (FDTD)*** 


*M.A. Fernandez and N. Zemzemi, 'Decoupled Time-Marching Schemes in Computational Cardiac Electrophysiology and ECG Numerical Simulation’, Mathematical Biosciences, Vol.226, No.|, pp.58-75, 2010. 
**F B. Sachse, C.D. Werner, K. Meyer-Waarden and O. Doessel, 'Solving the Forward Problem in Electrocardiography: Analysing the Selection of Appropriate Regions of Impedance Models’, Computers in 


Cardiology, Vol.25, pp.701-704, 1998. 
xM A. Stuchlyand T.W. Dawson, ‘Interaction of Low-Frequency Electric and Magnetic Fields with the Human Body’. Proceedings of the IEEE, Vol. 88, No. 5, pp. 643-664, May 2000. 
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2. VOLUME CONDUCTOR MODELING 
Finite Difference Method (FDM) 


Ease of implementation; 


Computational efficient; 


2D æ 3D =æ inhomogeneity 


Gradual model refinement. 


Surface potential 






+ |V dipole, lossless dielectric 
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+ |V dipole, lossless dielectric 


N 


2.VOLUME CONDUCTOR MODELING 


* US National Library of Medicine (NLM) dataset of cross-sectional 
photographs of human body 





e Correct dielectric characteristics assigned to all organs* 
*  Voxels (volume pixels) of Imm x Imm x Imm created 
* Visual Man body (parts) may be included in 3D FDM code 


*|. W. Massey, Creating AustinMan: An Electromagnetic Voxel Model of the Visible Human, Undergraduate Thesis, University of Texas at Austin, May 2011. 
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2. VOLUME CONDUCTOR MODELING 


Visible Human Project - AustinMan 


AustinMan is a voxel model of the human body 
that is being developed for electromagnetic 
simulations from the National Library of 
Medicine's Visual Human Man data set. The 
model is an “open model” that is easily 
accessible and verifiable. The model is freely 
available for research, teaching, or other non- 
commercial uses. 









AustinMan Electromagnetic Voxel models available at http://tiny.cc/AustinMan 
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2.VOLUME CONDUCTOR MODELING 
CST Voxel Family* 


CST Voxel Family 


c» сүт Anthropomorphic 
С һап models in 
resolutions ranging 
from 0.98x0.98x | Omm? 
to 2.08x2.08x8.0mm? 





| КАТЈА 
(pregnant) 


int 
EI 
[s 





www.cst.com | May-10 


*CST [Online]. Available at: http://www.cst.com/Content/Events/downloads/eugm20|0/Bio-Medical RF simulations with CST MICROWAVE STUDIO.pdf 
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Body surface electric field 


2. VOLUME CONDUCTOR MODELING 





Heart activity modeled as a non- 


rotating 60Hz sinusoidally excited 
dipole 


ИТС 





© IMEC 2015 


Torso modeled as consisting of heart, lungs bones and remaining tissue* 
*M.A. Fernandez and N. Zemzemi, ‘Decoupled Time-Marching Schemes in Computational Cardiac Electrophysiology and ECG Numerical Simulation’, Mathematical Biosciences, Vol.226, No.1, pp.58-75, 2010. 


3. EQUIVALENT VOLUME SOURCE MODEL 


e First body surface electric field results look promising; 
e Influence of body details don’t seem to be dominant; 


e Equivalent volume source model needed. 





Basic source volume models 





(Fixed) Dipole: single dipole, variable orientation and magnitude; 


e . Moving Dipole: as above, but also varying position; 
e Multiple Dipoles; 
e  Multipoles. 


Easy to implement; 
Accounts for 8076 of body surface potentials. 
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4.FIXED DIPOLE MODEL AND IMPLEMENTATION 
Equivalent Source Model 


e From ECG potential differences, measured over the Einthoven 
triangle leads, a source is constructed; 


e We try to find out if we can use this source in a CST voxel model 
to obtain ECG potentia differences. 





ATRIAL 





0 0 0 
80 1 0 
220 1 90 
230 2.5 0 
240 4 -30 
250 1:5 -80 
350 0 0 
450 2 15 
600 0 0 
1000 0 0 





' 
Sanada’ 


J ы 7 {А 
*J|. Malmivuo, R. Plonsey, Bioelectromagnetism, Principles and Applications of Bioelectric and Biomagnetic Fields, Oxford University Press, New York, 1995. 
ш 
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4.FIXED DIPOLE MODEL AND IMPLEMENTATION 


Equivalent Source Model 





- - Implementation in CST voxel model 

1 0 . 

| " e Voltage or current paths, having a 

2.5 0 е А . . 

à 3o sinusoidal excitation needed for 

"D performing a quasi-magnetostatic analysis; 
2 15 

2 о * Implement low-frequency sinusoidal 

0 0 


current path (dipole) with required 
amplitude and orientation per timestamp; 


e = Record maximum body surface electric field distribution; 


e Take differences for relative ECG signal strength. 


4 


Determine cross sectional heart plane in Gustav voxel model 
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4. FIXED DIPOLE MODEL AND IMPLEMENTATION 


Equivalent Dipole Source Plane and Orientations 


(95,45,405) 


(80,-50,405) 


(95,-65,405) 
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4. FIXED DIPOLE MODEL AND IMPLEMENTATION 


Equivalent Dipole Orientations 






dipole 
orientation 





8Oms 220ms 230ms 





240ms 250ms 350ms 450ms 600ms 





5. FIRST SIMULATION RESULTS 








5. FIRST SIMULATION RESULTS 
ECG Signal 





Difference in 60Hz maximum body surface potential values 


Fixed dipole position (sinus node in heart), varying amplitude and direction 


Б 


PQRST complex 





Relative amplitude 








0 100 200 300 400 500 600 
Time (ms) 
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6. CONCLUSION 


To use a CST human voxel model for simulating ECG-like signals on 
the body 


e A fixed dipole equivalent volume source is needed; 
* This source is implemented as a sinusoidally excited current path; 


* The maximum body surface electric field is recorded for different time 
stamps; 


• Тһе current path is rotated and changed in magnitude for every time stamp. 
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